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Josephson junctions containing muhiple ferromagnetic layers can carry spin-triplet supercurrent 
under certain conditions. Large-area junctions containing multiple domains are expected to have 
a random distribution of or tt coupling across the junction surface, whereas magnetized samples 
should have uniquely tt coupling everywhere. We have measured the area dependence of the critical 
current in such junctions, and confirm that the critical current scales linearly with area in magnetized 
junctions. For as-grown (multi-domain) samples, the results are mixed. Samples grown on a thick 
Nb base exhibit critical currents that scale sub-linearly with area, while samples grown on a smoother 
Nb/Al multilayer base exhibit critical currents that scale linearly with area. The latter results are 
consistent with a theoretical picture due to Zyuzin and Spivak that predicts that the as-grown 
samples should have global 7r/2 coupling. 

PACS numbers: 74.50.-|-r, 74.45.-l-c, 75.70.Cn, 74.20.Rp 



I. INTRODUCTION 

When a superconducting (S) metal is placed in contact 
with a nonsuperconducting (N=normal) metal. Cooper 
pairs can "leak" out of the superconductor and modify 
the properties of both materials. This process, called the 
superconducting proximity effect, has been studied for 
several decadesji When the normal metal is replaced by 
a ferromagnetic (F) metal, the two electrons of the pair 
enter different spin bands with different Fermi wavevec- 
tors. As a result, the pair correlations oscillate and decay 
rapidly with increasing distance from the S/F interfacCi^ 
The oscillating, short-ranged proximity effect in S /F sys- 
tems was predicted as early as 1982;^ and has been ob- 
served convincing by many groups over the past decade 

In 2001, it was predicted that pair correlations with 
spin-triplet symmetry can be generated in S/F systems 
containing certain kinds of magnetic inhomogeneities in- 
volving non-collinear magnetizations, even if all of the 
superconducting materials in the system have conven- 
tional spin-singlet symmetry..Sr— Spin-triplet pair corre- 
lations do not experience the exchange field in F, hence 
the proximity effect due to those pair correlations is long 
ranged. Some experimental evidence for spin-triplet cor- 
relations appeared in 2006^^^ then more convincing ev- 
idence appeared in 2010iii^— Our own contributioi>ii 
was based on measurements of the critical current Ic in 
Josephson junctions of the form S/F7SAF/F"/S, where 
SAF stands for "synthetic antiferromagnet" and F' and 
F" are thin ferromagnetic layers whose magnetizations 
must be at least partly non-collinear with that of the 
SAF. The SAF is a Co/Ru/Co trilayer with Ru thick- 
ness 0.6 nm, which causes anti-parallel coupling of the 
two surrounding Co magnetizations. Because of that 
anti-parallel coupling, the SAF produces nearly zero net 
magnetic flux, which would otherwise distort the "Fraun- 
hofer patterns" one observes in plots of Ic vs magnetic 
field H applied in the plane of the Josephson junction. 
We found that Ic in our samples hardly decreased as the 
total Co thickness was increased up to 30 nm, whereas 



Ic decayed rapidly in similar samples without the F' and 
F" layersji^ The long-range nature of the supercurrent in 
the samples containing F' and F" layers provided strong 
evidence for the spin-triplet symmetry of the current- 
carrying electron pairs. Furthermore, we have recently 
shown that Ic increases further when we magnetize our 
samples with an in-plane applied magnetic fieldi^ The 
explanation is that the F' and F" layers are magnetized 
parallel to the field, while the SAF undergoes a "spin- 
flop" transition whereby the two Co layers end up with 
their magnetization perpendicular to the direction of the 
applied field.—— According to theory,—"— this configu- 
ration with perpendicular magnetizations maximizes the 
magnitude of the spin-triplet supercurrent. 




FIG. 1: Cartoon showing relative orientations of magnetiza- 
tion for the ferromagnetic layers in our Josephson junctions. 
If angles 8i and 62 have the same sign (where we constrain 
^i|, 1^2! < tt), the junction will have tt coupling; if they have 
opposite signs, the junction will have coupling. 

The results described above raise several questions, 
only some of which have been answered by subsequent 
work in our groupj^^^ The question that motivated this 
paper arises from the theoretical prediction that Joseph- 
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son junctions of the form S/F'/F/F"/S, carrying spin- 
triplet supercurrent, can be either in the 0-state or the 
TT-state depending on the relative orientations of the three 
ferromagnetic Iayersfi2i2£ (It does not matter whether 
the central F layer is a single ferromagnetic layer or an 
SAFj2ii2^) The situation is illustrated in Figure [U which 
shows the relative orientations of the magnetizations of 
all four ferromagnetic layers in our junctions. According 
to theory, if the two angles di and 62 have the same sign, 
then the junction will have tt coupling; if they have op- 
posite signs, the junction will have couplingi^"— (We 
define the angles by the constraint 16*11, 16*21 < tt.) Since 
the magnetic layers in our samples consist of many do- 
mains when the samples are first grown, we would ex- 
pect the Josephson coupling in our junctions to exhibit 
a random spatially- varying pattern of 0-coupling and tt- 
coupling across the junction area. In that case, if a fixed 
difference in gauge-invariant phase is applied across the 
junction, some areas of the junction will provide positive 
supercurrent while others will provide negative supercur- 
rent - i.e. supercurrent flowing in the opposite direction. 
One could then calculate the total supercurrent naively 
using an analogy to the random walk problem: while the 
mean supercurrent averaged over many domains would 
be zero, the typical supercurrent in a given sample would 
be proportional to the square-root of the number of do- 
mains, hence to the square-root of the junction area. Af- 
ter the samples are magnetized, the magnetizations of 
the F' and F" layers are parallel to each other, hence 6*1 
and 62 have the same sign and the junction should have 
TT-coupling everywhere. In that case, the critical super- 
current will be proportional to the junction area, as is 
the case in conventional Josephson junctions. 

A completely different view of a Josephson junction 
containing a random spatially-varying pattern of and 
TT couplings has been proposed by Zyuzin and Spivak 
(ZS)i^ Those authors addressed S/F/S junctions with 
spin-singlet rather than spin-triplet supercurrent, and 
considered the situation where the F-layer thickness is 
large, so that the average supercurrent is small, whereas 
mesoscopic fluctuations of the Josephson coupling have 
random sign. In our samples, the spin-singlet supercur- 
rent is negligibly small (c.f. Figure 3 in Ref. and 
the random-sign spin-triplet Josephson coupling arises 
from the local variations in magnetic domain struc- 
ture. In spite of the different mechanisms underlying the 
spatially-varying random-sign Josephson coupling, there 
is no apparent reason why the ZS model should not ap- 
ply to our spin-triplet Josephson junctions. ZS calcu- 
lated the total energy of such a junction, and concluded 
that the ground state corresponds to, on average, a 7r/2 
phase difference between the two superconducting elec- 
trodes. The phase difference is spatially modulated, with 
local variations toward lower phase difference in regions 
of 0-coupling and larger phase difference in regions of 
TT-coupling. According to the ZS result, the total super- 
current scales with the junction area, as is the case for 
conventional Josephson junctions. 



The purpose of this paper is to measure experimentally 
the area dependence of the supercurrent in our Josephson 
junctions, to determine which of the pictures presented 
above applies. 



II. SAMPLE FABRICATION 
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FIG. 2: (color online) Schematic diagram of Josephson junc- 
tion samples (not to scale). The current flows in the vertical 
direction. 

The Josephson junctions in this work were fabricated 
by dc triode sputtering, photolithography and ion- 
milling, as described in our previous publications J^ii^ 
The structure of the junctions is shown schematically 
in Figure [2j In this work, we have grown two types 
of samples with different superconducting base lay- 
ers: either a single 150-nm layer of Nb, or a Nb/Al 
multilayer stack described below. For the first kind of 
sample we start by growing a multilayer of the form 
Nb(150)/Cu(5)/Ni(1.5)/Cu(10/Co(6)/Ru(0.6)/Co(6)/ 
Cu(10)/Ni(1.5)/Cu(5)/Nb(20)/Au(15), where all thick- 
nesses are in nm. That stack is sputtered in one run 
without breaking vacuum. For the second type of sample, 
the Nb base layer was replaced with a Nb/Al multilayer of 
the form [Nb(40nm)/Al(2.4nm)]3/Nb(40nm)/Au(15nm). 
This structure was motivated by the work of Thomas et 
al»^, who used a thin Nb/Al multilayer on top of thick 
Nb to reduce surface roughness. For those samples, the 
chamber was opened briefly after the Au deposition, 
while flowing N2 gas, to change sputtering targets. (Our 
sputtering system holds 6 targets, whereas the second 
type of samples require 7 different materials.) After 
pumpdown, we sputtered first 20 nm of Nb, followed by 
the rest of the stack up through the top Au layer. We 
know from our top electrode fabrication procedure that 
the top 15-nm Au layer adequately protects the under- 
lying Nb from oxidation, and is driven superconducting 
when sandwiched between two Nb layers. The same 
should be true for the bottom Au layer protecting the 
Nb/Al multilayer base. 
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III. EXPERIMENTAL RESULTS 




FIG. 3: (color online) Atomic force microscopy pictures of (a) 
a 200- nm thick Nb base layer and (b) a Nb/Al multilayer as 
described in the text. 



To ascertain the surface roughness of the two types 
of base layers, we performed atomic force microscopy 
(AFM) measurements on a bare 200-nm Nb base and 
on a Nb/Al multilayer stack (up to the Au layer dis- 
cussed above). The results are shown in Figure [3l The 
root-mean-squared roughnesses of the first and second 
base layers are 0.53 nm and 0.23 nm, respectively, over 
the 250 X 250nm^ area shown. As expected, the Nb/Al 
multilayer provides a smoother base than the pure Nb. 



All of the data reported here were acquired at 4.2K 
with the sample dipped into a liquid helium dewar. 
A current comparitor circuit using a Superconducting 
Quantum Interference Device (SQUID) as a null detector 
was used to measure the current-voltage (I-V) character- 
istic of the samples.— All samples exhibit the standard 
I-V characteristic for an overdamped Josephson junction: 



Vil) = Sign[I] * i?ArRe[(/2 - /2)i/2] 



(1) 



where i?jv is the normal-state resistance determined from 
the slope of the V-I relation at large currents. 



Fraunhofer patterns and result of magnetizing 
samples 




FIG. 4: (color online) Critical current vs. applied field for 
two 3-fim diameter Josephson junctions, measured in the vir- 
gin state (panels a and c) , and after the samples were magne- 
tized by a large in- plane field (panels b and d). In the virgin 
state, two separate runs are shown for each sample. The lines 
connect the data points; they are only guides for the eye. 



For both types of samples, Josephson junctions with 
circular cross section were defined by photolithography 
and ion milling. Because we wanted to obtain data on 
multi-domain samples covering a large dynamic range of 
areas, we fabricated junctions with diameters of 3, 6, 12, 
24, and 48 micrometers. Unfortunately, the largest sam- 
ples rarely produced high-quality data, hence we restrict 
ourselves here to the samples of diameters 3, 6, and 12 
fim. Another difference between the samples measured 
here and those measured in our previous publications is 
that these were ion milled only to the top copper layer 
in order to keep the magnetic domain structure intact, 
whereas our previous samples were typically ion milled 
partway through the top Co layer. 



All samples are initially characterized by applying a 
small magnetic field perpendicular to the current direc- 
tion, i.e. in the plane of the substrate. A plot of Ic 
vs H should yield the classic "Fraunhofer pattern" (ac- 
tually an Airy pattern for our circular pillars). Figure 2] 
shows representative Ic vs H data for two 3-/im-diameter 
junctions, in the virgin state (panels a and c), and after 
being magnetized by a large in-plane field (panels b and 
d). We measured a few samples several times to deter- 
mine how much the data vary from run to run; panels a 
and c show two virgin-state runs for two of these sam- 
ples. Several features are evident from the data: i) the 
Frauhofer patterns in the virgin state fluctuate from run 
to run; ii) the quality of the Fraunhofer patterns is better 
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after the samples are magnetized than in the virgin state; 
iii) Ic is enhanced after the samples are magnetized, as 
we reported recentlyf^^ and iv) the central peak in the 
Fraunhofer patterns of the magnetized samples is shifted 
to negative field by about 30 Oe. 

The variability and relatively low quality of the 
Frauhofer patterns in the virgin state are undoubtedly 
due to the random domain structures of the ferromag- 
netic layers in the samples.— One cannot know a pri- 
ori if the problem is due to the 1.5-nm thick Ni F' and 
F" layers or to the two 6-nm thick Co layers making up 
the SAF. We believe it is the former, given the improve- 
ment in the quality and reproducibility of the Fraunhofer 
patterns after magnetization. Magnetizing the samples 
forces the Ni domain magnetizations to point in nearly 
the same direction, and probably causes the average do- 
main size to grow. In contrast, the Co/Ru/Co SAF is 
expected to have its best antiparallel coupling in the vir- 
gin state. 

The enhancement of Ic by magnetizing the samples has 
potential contributions from two factors. In our recent 
work,'^^ we emphasized optimization of the angles 6i and 
O2 between the inner and outer ferromagnetic layer mag- 
netizations. Those angles vary randomly across the junc- 
tion area in the virgin-state samples, whereas they should 
both be close to the optimal value of 7r/2 after the sam- 
ples are magnetized, due to the Co/Ru/Co SAF under- 
going a spin-flop transition. A second contributor may 
be the fact that, in the virgin state, the Josephson cou- 
pling varies randomly between 0-coupling and 7r-coupling 
across the junction area. If the random- walk picture dis- 
cussed earlier is valid, then one would expect the value 
of Ic in a typical sample to scale with the square-root of 
the junction area, as discussed earlier. After the samples 
are magnetized, there should be 7r-coupling everywhere, 
so that the supercurrent adds constructively across the 
entire junction area. 

Finally, the shift in the Fraunhofer patterns after mag- 
netization has been observed previously;^"— The central 
peak occurs at the point where the flux due to the ap- 
plied field exactly cancels the intrinsic flux due to the 
magnetization of the Ni layers inside the junction. 

The evolution of Ic as the sample is magnetized is 
shown for a 6-/im diameter Josephson junction in Figure 
[5l The sample was first measured in the as-grown state 
(H = 0). Then the magnetizing field H was stepped up 
to 3600 Oe with varying step sizes evident in the figure. 
After application of each value of the field is reduced 
to zero and the Fraunhofer pattern is measured in low 
field. The squares show the resulting values of IcRn as 
the sample is magnetized. For low fields, nothing hap- 
pens. Then there is a shallow dip in IcRn for H near 500 
Oe. That dip is observed in many samples, but is not 
fully understood; it may be related to a change of the Ni 
domain structure. When H is increased above 500 Oe, 
Ic increases sharply. The field range where Ic increases 
corresponds to the field range where the Ni films become 
magnetized. (We know this because the field where Ic 
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FIG. 5: (color online) IcRn product vs. applied in-plane field 
for a 6-/ini diameter Josephson junction. The sample was 
first magnetized in positive field (squares), then the sample 
was demagnetized and finally re-magnetized in negative field 
(circles). 

increases varies with Ni layer thickness, and matches the 
coercive field of the Ni determined from separate magne- 
tization measurements of large-area Ni filmsJ^) Figure 
[5] also shows what happens when a field is applied in 
the opposite direction to the original magnetizing field 
(circles). Again, nothing happens for small field values. 
Then, as the Ni films are demagnetized, Ic drops to val- 
ues as low as or even lower than the value at the dip we 
observed when first magnetizing the samples. As the Ni 
films are re-magnetized in the negative direction, Ic in- 
creases sharply again to a value essentially identical with 
that observed on the positive field side. We have mea- 
sured full magnetization curves for several samples, and 
they all look very similar to the one shown in Figure [5j 



B. Area dependence: samples with Nb base layer 

To provide good statistics and to reveal the extent of 
sample-to-sample fluctuations, we have fabricated and 
measured a large number of Josephson junctions with 
diameters of 3, 6, and 12 /im. Figure [5] shows the re- 
sults for samples grown on our traditional thick (150 
nm) superconducting Nb base, for both the virgin and 
magnetized states. The points representing the magne- 
tized state (open symbols) are averages of the measure- 
ments taken after application of 1600, 2000, and 2400 
Oe. (There is little variation of IcRn between those three 
measurements.) The points representing the virgin state 
(solid symbols) are usually averaged over two runs, al- 
though a few samples were measured only once, and one 
sample was measured 5 times in the virgin state. 

The results of the magnetized state measurements in 
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FIG. 6: (color online) Critical current times normal-state re- 
sistance vs. junction diameter for Josephson junctions grown 
on a 150-nm Nb base electrode. Solid symbols represent 
virgin-state data; open symbols represent data acquired after 
the samples were magnetized by a large (~2000 Oe) in-plane 
magnetic field. 



Figure [S] clearly show that IcRn is essentially indepen- 
dent of sample area. Since Rn is inversely proportional 
to junction area, this means that Ic is proportional to 
area. That is the usual situation, and is what one ex- 
pects when the Josephson coupling is uniform across the 
junction area. In contrast, the virgin-state data show a 
decrease in IcRn with increasing sample size. According 
to the random walk model discussed in the introduction, 
Ic should scale with the square-root of the junction area, 
hence IcRn should scale inversely with the square-root 
of area, or equivalently, inversely with junction diameter 
D. The virgin-state data shown in Figure [5] do exhibit 
a noticeable decrease with junction diameter, support- 
ing the random walk picture, although the dependence is 
slightly less steep than Ic oc D^^ . 



C. Samples with Nb/Al base layer 

As shown in Figure [31 the Nb/Al multilayer base pro- 
vides a smoother surface than the pure Nb base layer. We 
were curious as to whether the smoother base would influ- 
ence any of the Josephson junctions properties. We per- 
formed the same measurements of Ic on the second batch 
of samples, grown on the smoother Nb/Al base, as were 
performed on the first batch, grown on pure Nb. The re- 
sults are shown in Figure [T] In the magnetized state, the 
results agree closely with those of the first batch, shown 
in Figure ini Not only is IcRn essentially independent 
of junction diameter, but the actual values of IcRn are 
very close to those of the previous batch of samples. Af- 
ter being magnetized, there is remarkable consistently in 
the values of IcRn for the 21 samples displayed in Fig- 
ures [B] and [71 In the virgin state, however, the story is 
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FIG. 7: (color online) Critical current times normal-state re- 
sistance vs. junction diameter for Josephson junctions grown 
on a Nb/Al multilayer as described in the text. Solid sym- 
bols represent virgin-state data; open symbols represent data 
acquired after the samples were magnetized by a large ( 2000 
Oe) in-plane magnetic field. 

different. In contrast to what we observed in the first 
batch of samples, IcRn in the second batch hardly varies 
with junction diameter. These samples deposited on top 
of the smoother base electrode appear to validate the ZS 
theory, which predicts that Ic should scale linearly with 
junction area. 




FIG. 8: (color online) Summary of IcRn data for all the 
Josephson junctions studied in this work. Each symbol rep- 
resents the average value for all samples of a given size and 
base layer, in either the virgin-state (solid symbols) or af- 
ter being magnetized (open symbols). The circles represent 
samples grown on a 150-nm Nb base layer, while the tri- 
angles represent samples grown on a Nb/Al multilayer de- 
scribed in the text. The dot-dashed line illustrates the rela- 
tion IcRn oc . 

The situation is summarized in Figure [HI where we 
have averaged together the values of IcRn for all sam- 
ples of a given diameter, fabricated on a given base layer. 
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With some of the sample-to-sample fluctuations averaged 
out, the trends are striking: i) The magnetized data are 
remarkably consistent, and hardly depend on the base 
layer; ii) the virgin-state data from samples deposited on 
the thick Nb base exhibit IcRn values that decrease sub- 
stantially with junction diameter, but not quite as fast 
as , which is shown by the dot-dashed line in the fig- 
ure; iii) the virgin-state data from samples deposited on 
the smoother Nb/Al multilayer base exhibit IcRn values 
that are independent of junction diameter. 

What are we to make of the results shown in Figure |5]? 
It appears that the roughness of the base layer has a pro- 
found effect on the area scaling of Ic in the virgin state. 
One can imagine several possible explanations. One is 
that the spectrum of Andreev boimd states, and hence 
the Josephson current, in an S/F/S junction depends in a 
fundamental way on whether the S /F interface is smooth 
or rough. There has been some theoretical work on inter- 
face roughness in S/N systems, "^^ but we are not aware 
of any such work that addresses our results directly. A 
second possibility is that the roughness of the Nb base 
layer perturbs the domain structure of the Ni F' and F" 
layers - possibly even to the extent that one or both of 
those layers are not continuous. If this is the case, one 
still has to explain how modified or discontinuous F' and 
F" layers would affect the area scaling of the junction 
critical current. We know that the coercive fields of the 
Ni layers increase continuously with decreasing Ni thick- 
ness down to 1 nm, as shown in the data in Figure 2 of 
Ref. 16], so there does not appear to be any anoma- 
lous magnetic behavior induced by the rough base. A 
third possibility is that the observed sub-linear scaling 
of the supercurrent with junction area for the junctions 
grown on the rougher Nb base is simply a reflection of 
the gradual deterioration of the quality of the Fraunhofcr 
patterns with increasing sample size. A possible way to 
ameliorate that issue in the future would be to use PdNi 
alloy rather than pure Ni as the F' and F" layers. PdNi 
is a weak ferromagnetic material with small magnetiza- 
tion, and in earlier work we were able to produce Joseph- 
son junctions with high-quality Fraunhofer patterns even 
with much thicker PdNi layers than one would need for 
this experiment f2§. The optimal PdNi thickness for pro- 



ducing spin-triplet supercurrent is in the range of 4-6 nm, 
which is much thicker than the 1-2 nm optimal range for 
Niiiii2^ We chose pure Ni for the F' and F" layers in the 
present work because it produces the largest values of Ic 
in the virgin state,— but thicker PdNi layers might be 
less sensitive to the nm-scale roughness of the Nb base. 



IV. CONCLUSIONS AND OUTLOOK 

In summary, we have measured the area-dependence 
of the critical current in S/F/S Josephson junctions car- 
rying spin-triplet supercurrent. After the samples are 
magnetized, the critical current has its largest value, and 
it scales linearly with area as is the case for conventional 
Josephson junctions. In the virgin state, however, the 
results are mixed. Samples grown on our traditional 
thick Nb base exhibit critical currents that grow sub- 
linearly with area, whereas samples grown on a smoother 
Nb/Al multilayer base exhibit critical currents with con- 
ventional area scaling. The former may be an indication 
that the supercurrent is not uniform over the junction 
area, while the latter provides indirect support for a the- 
oretical model of Zyuzin and Spivak 

In the future, it would be interesting to measure the 
current-phase relation of our junctions, both in the vir- 
gin and magnetized states. According to the ZS theory;^ 
the virgin state junctions should be in a 7r/2 state, while 
according to spin-triplet junction theory^^^— the magne- 
tized junctions should be in the tt state. Current-phase 
measurements are technically more challenging than the 
critical current measurements reported here, but they 
might provide more direct evidence of the underlying 
physics than do the area-dependent measurements re- 
ported area. 
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